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Quantum tunneling in an asymmetric (with strongly different capacitances) SQUID is studied. 
Since capacitances play a role of masses one phase, related to a large mass, becomes "heavy" and 
remains always a constant in a tunneling process. Tunneling in an asymmetric SQUID becomes 
one-dimensional with a condition of optimization of tunneling probability with respect to a value of 
the "heavy" phase. An unusual temperature dependence of the tunneling probability is obtained. 
It has a finite slope at zero temperature and a transition between thermally assisted tunneling and 
pure activation can be not smooth depending on current through a SQUID. 
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I. INTRODUCTION 

A decay of a zero voltage state via quantum tunnel- 
ing of a phase across a potential barrier in Josephson 
junctions is possiblcii"— Tunneling in a single Joseph- 
son junction is similar to a conventional one-dimensional 
quantum mechanical process. In this case the tunneling 
mechanism is described by theory of Wentzel, Kramers, 
and Brillouin (WKB).^ Tunneling occurs from a classi- 
cally allowed region which is a conventional potential well 
where energy levels are quantized.—^— Quantum coherence 
between potential wells was demonstrated J^^J^ 

Besides single Josephson junctions, superconducting 
quantum interference devices (SQUID) are also a mat- 
ter of active investigation.—"— A SQUID consists of two 
junctions and, therefore, represents a two-dimensional 
system where macroscopic quantum tunneling is also pos- 
sible. 

Tunneling in multi-dimensional systems was described 
in literature.—"— According to those theories, there is a 
certain underbarrier path where a wave function is local- 
ized and it decays along the path. However, a tunneling 
scenario in two (many) dimensions can substantially dif- 
fer from that main-path mechanism. For example, in a 
SQUID a main path can split by two onesi^ and even 
by an infinite number of equivalent paths which interfer 
providing multi-path tunneling. See also2§r— . 

Capacitances of two junctions in a SQUID play a role 
of masses. When one capacitance is large the masses 
become very different and one phase becomes "heavy". 
This provides an additional interest for study of an asym- 
metric SQUID. A behavior of an asymmetric SQUID af- 
ter tunneling was investigatedi^ In this paper we also 
focus on strongly asymmetric SQUIDs, namely, on tun- 
neling process in them. 

During tunneling process a motion along the "heavy" 
coordinate is weakly generated and the process of bar- 
rier crossing becomes almost one-dimensional when the 
"heavy" phase is fixed. This fixed value should be de- 
termined from a condition of maximum of a tunneling 
probability. 



That program has been performed in the paper. There 
are two unusual features of results. 

First, the tunneling probability, as a function of tem- 
perature, has a finite slope at low temperature. This con- 
trasts to a temperature dependence for a one-dimensional 
barrier where that slope is zero. 

Second, a transition at a finite temperature between 
thermally assisted tunneling and pure activation changes 
its character when current approaches the critical value. 
At those currents temperature dependence of tunneling 
probability exhibits a finite jump of slopes at the transi- 
tion temperature. When current is not too close to the 
critical value the transition is smooth as for a one dimen- 
sional barrier. 

In Sec. mil we apply to a SQUID a semiclassical for- 
malism of Hamilton- Jacobi. In Sec. IIVI the method of 
classical trajectories in imaginary time is used which ac- 
counts for an optimization of tunneling probability with 
respect to a value of the "heavy" phase. In Sec. IVIII it 
is argued that experimental observations of the proposed 
phenomena in SQUID is real. 



II. FORMULATION OF THE PROBLEM 

We consider a dc SQUID, consisting of two Joseph- 
son junctions with phases (pi and ip2, with no dissipation 
when the two junctions are inductively coupled. Critical 
currents of the junctions are equal but capacitances, Ci 
and C2, are different so that 
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FIG. 1: (Color online) Curves of a constant energy V{x, y) = 
E ai a — 0.9. Tunneling occurs along the dashed line, y ~ yo, 
where an underbarrier wave function is localized. 



where the dimensionless current j ~ I/'^Ia the Joseph- 
son energy Ej — hlc/2e, the plasma frequency uj = 
y/2eIc/hCi, and the coupling parameter /3 = 2'kLIc/^o 
are introduced. Here Ic and L are critical current and in- 
ductance of each individual junction. The magnetic flux 
quantum is $o = Trfic/e. As follows from Eqs. ([T]) and 
(12), (^2 is a "heavy" phase. 

Below we consider a large /3 and the total current / 
close to its critical value, (1 — j) <C 1. New variables are 
introduced by the relations 



^1 = 1 + (3^-i)v/2(r^+^ 

¥^2 = 1 + (3^-1)^2(1^+^. 



Below time is measured in the unit of 

u) \ 1 + a 
where the coupling parameter is 
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The energy 



takes the form 
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The potential energy is 

V{x,y) = Vo{x) + Vo{y) 



+ Vix,y) 

2axy 
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where Vo{x) = x^ — x^. B in Eq. ^ is called semiclas- 
sical parameter. When B is large the phase dynamics is 
mainly classical. Below we consider that case, 1 <^ B. 
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FIG. 2: (a) Form of 5E{y). (b) Corresponding density distri- 
bution is Gaussian. It is plotted for some x under the barrier. 



III. DESCRIPTION OF TUNNELING 

A classical dynamics of phases in a SQUID relates to 
Eqs. ^ and ([H]). The effective particle moves in the clas- 
sically allowed region, in a vicinity of the point x = y = 0, 
which is restricted by the potential barrier. As known, 
the particle can tunnel through the barrier resulting in 
experimentally observable phase jumps. Character of 
tunneling depends on coupling strength a between the 
two junctions. At a = 0.90 the curves of equal potential, 
V{x, y) ~ E, are shown in Fig. [TJ An effective particle 
tunnels from one classically allowed region (the potential 
well) to another (the outer region). 

To quantitatively study the problem of two- 
dimensional tunneling one should solve the Schrodinger 
equation with the exact potential (|S]). Since the potential 
barrier is almost classical one can apply a semiclassical 
method when a wave function has the form 
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where the classical action is hBa. a satisfies the equation 
of Hamilton- Jacobi^ 
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We define the energy E by the relation 
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At a large M a solution of Eq. ((TU)) can be written in the 
form a = ao + ai where ai is small and (Tq is given by 
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where 5E{x, y) is some function to be specified. It is easy 
to conclude that the correction cri is small (proportional 
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FIG. 3: (a) Form of SE{y). (b) Corresponding density distri- 
bution, plotted for some x under the barrier, is not Gaussian. 
See the text. 
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FIG. 4: (Color online) F is tunneling probability and tem- 
perature T is measured in the units of hu/y/p. The numbers 
mark values of the parameter (1 — j)l3'^ . 



probability of tunneling with a fixed energy E is 



to 1/y/M) when the derivative d5E/dx is small (propor- 
tional to So we consider below 5E{y). 

The function 5E(y) is determined by a state in the well 
from which tunneling occurs. 

When the case of Fig. is realized the last term 
in Eq. provides a Gaussian distribution of density 
around the line y — yo shown in Fig. ^ih) for some x 
under the barrier. This is analogous to a conventional 
scenario of tunneling in two dimensions >2^i2S 

In the case of 8E{y) of Fig. [3ja) at yo < y the density 
drops down under the barrier due to the second term in 
Eq. (fT2|) and at y < yo due to the first one. Therefore 
the underbarrier density in that case is also localized in 
a vicinity of the line y = yo- It is shown in Fig. |3ljb) 
for some x under the barrier. This is not a Gaussian 
distribution but one of the type exp[— c(j/ — yo)^^'^]- At 
yo < y the parameter c ^ %/M is large. 

The situations in Figs. [2] and [3] relate to different types 
of states in the potential well from which tunneling oc- 
curs. In the case of Fig. [2] a Gaussian distribution of 
density holds also in the well because the last term in 
Eq. (fT2|) dominates. In the case of Fig. [3] the distribu- 
tion in the well is analogous to Fig. EJb) when the part 
at y < yo is horizontal and, therefore, the state is dis- 
tributed over a finite distance y in the well. 

In the both cases. Figs. [5] and 121 tunneling occurs along 
the certain line y = yo and yo should be determined 
from the condition of maximum of a tunneling proba- 
bility. This corresponds to classical mechanics when a 
particle does not move along a "heavy" direction. With 
an exponential accuracy a tunneling probability is the 
same for the both types of states in the well. 



IV. TUNNELING PROBABILITY 

Since tunneling occurs along the line y = yo one can 
use a WKB approach as in a one-dimensional case. The 
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where 



The one-dimensional potential v{x) is given by 



v{x) = X — X 
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The integration in Eq. (fT4|) is restricted by the classically 
forbidden region where E < v{x). 

Tunneling probability at a fixed temperature accounts 
for the Gibbs factor and is determined by 
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with a subsequent optimization with respect to E. Tak- 
ing Eq. (|lip , one can write Eq. (|16p in the form 

r ~ exp {^2BAt) , (17) 

where 

At^A{E) + j. (18) 
The parameter 9 is connected with temperature 
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Minimization of At with respect to energy defines E 
by the equation 
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y/2j ^v{x) -E' 
The parameter yo should be chosen to minimize At- By 
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FIG. 5: (Color online) Amplification of the lower set of curves 
in Fig.|4] Left parts of the curves relate to a thermally assisted 
tunneling and right parts pertain to pure activation. 



introducing imaginary time t ■ 
written in the form 



ir the action At can be 



At = 
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dT 



dx 
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where the classical trajectory under the barrier is deter- 
mined by Newton's equation 
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with zero velocities at the terminal points, r = and 
T = 1/6. According to Eq. (PO]) . 1/9 is the underbarrier 
time of motion between two terminal points. In terms of 
trajectories, the condition of minimum At with respect 
to yo takes the form 
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For a strongly asymmetric SQUID, a large M, tun- 
neling occurs along a straight line y = yo shown in 
Fig. [TJ The action I^T^ depends on two parameters, a 
and T^/P/hjj. 

The tunneling probability satisfies the relation 
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A recipe of calculation of the action At is the following. 
At fixed a, yo, and 9 one should find a solution of Eq. (|22|) 
with zero velocities, dx/dr = 0, at r = 0, 1/0. That 
solution has to be inserted into the relation p3| which 
defines j/o at fixed 9 and a. The solution with the defined 
?/o should be substituted into Eq. (PT|) which produces 
At {a,Ty/0/huj). We demonstrate in Sec. IVl how this 
scheme works in the case of low temperatures. 



V. TUNNELING AT LOW TEMPERATURES 

At low temperatures the energy E should be close to 
the minimum of the potential v{x) providing a long un- 
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the action takes the form 
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A minimization with respect to j/q of the action (|26p . ac- 
counting for ((25)) . is equivalent to Eq. ((23|) . The resulting 
action, at low dimensionless temperature 0, is 
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VI. RESULTS 



We performed a numerical solution of Eq. p2|) . The re- 
sults for the tunneling probability are presented in Figs.|3] 
and [5] where temperature T is measured in the units of 
hbj/ \f^. Each curve in Figs.|4]and[5]consists of two parts. 
To the left of a dot each curve relates to above trajectory 
calculations corresponding to thermally assisted tunnel- 
ing. To the right of a dot a curve is solely due to thermal 
activation 
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The activation energy is given by the saddle point 
Vixs^ys) which coinsides with the crossing point of two 
curves in Fig.[T] The steepest descent in Fig.[T]goes along 
the direction x — y. The saddle point {xs,ys} is deter- 
mined by the conditions dV{x, y) /dx = dV{x, y) /dy = 0. 
At 1/2 < a 
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At a < 1/2 the transition to a pure activation regime 
is smooth as in a one-dimensional case. It is analogous 
to type II phase transition. This corresponds to 2 < 
(1 - in Figs. S] and [S] At 1/2 < a the transition 
to the activation regime reminds type I phase transition. 
A derivative with respect to temperature jumps at those 
points which can be observed in Fig. m 

Numerically calculated curves in Figs. |4] and [5] match 
at low temperatures the analytical dependence followed 
from Eqs. (l24l) and ((27|) . At low temperatures 
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FIG. 6: (Color online) Tunneling probability F versus current 
for various values of the dimensionless temperature T^/p/hio. 
The dashed curve separates thermally assisted tunneling from 
a pure thermal activation. 



We note that the slope in the temperature dependence 
([5T|) is finite. 

The tunnehng probability F as a function of the pa- 
rameter (1 — is plotted in Fig. [6] for different values 
of the dimensionless temperature T\fp/huj. This plot 
shows how F depends on current at a fixed temperature. 

VII. DISCUSSIONS 

Quantum tunneling across a one-dimensional static po- 
tential barrier is described by WKB theory. Accord- 
ingly, in two dimensions the main contribution to tunel- 
ing probability comes from an extreme path linking two 
classically allowed regions. The path is a classical tra- 
jectory with real coordinates which can be parametrized 
by imaginary time. The under barrier trajectory is a so- 
lution of Newton's equation in imaginary time. Under 
the barrier the probability density reaches a maximum 
at each point of the trajectory along the orthogonal di- 
rection with respect to it. Along that direction the den- 
sity has a Gaussian distribution. Therefore around the 
trajectory, which plays a role of a saddle point, quantum 
fluctuations are weak. The wave function, tracked along 
that trajectory under the barrier, exhibits an exponen- 
tial decay analogous to WKB behavior. This constitutes 
a conventional scenario of tunneling in multi-dimensional 
cas o^^'^^ which can be called main-path tunneling.^^ 

In a symmetric (not very asymmetric) SQUID besides 
the conventional main-path tunneling^^ also multi-path 
tunneling is possiblci^ In that case a density distribution 
under a barrier is not as in Fig.[2][b), but of the type as in 
Fig. |3l[b). For a very asymmetric SQUID, considered in 
this paper, the both mechanisms result in the same tun- 
neling exponents since the underbarrier channel shrinks 
to the line y = yo due to the mass difference. 



We used a semiclassical approximation when there are 
many levels in the well. This approach sometimes is 
not appropriate in one-dimensional Josephson junctions 
where a barrier is weakly transparent but nevertheless 
there is only a few levels in the well (say, five) ^dn^i^ In a 
SQUID based on two such junctions the number of lev- 
els can be roughly estimated as 5 x 5. In our case of 
a strongly asymmetric SQUID that number should be 
multiplied by the large parameter \/M. Therefore the 
approximation of a large number of levels in the well is 
reasonable. 

We propose two peculiarities of tunneling in an asym- 
metric SQUID which do not exist in a single junction and 
in a not very asymmetric SQUID. 

One of them is temperature dependence of tunneling 
probability at low temperature. According to Eq. (j3ip . 
the curves in Fig.[5]have a finite slope at low temperature. 
In one dimension the slope is zero due to the exponent 
exp{-const/T) instead of T in Eq. (|3T|) . 

The second peculiarity is an unusual transition be- 
tween thermally assisted tunneling and pure activation 
marked by dots in Fig.[5l At 2 < (1 — j)/3^ the transition 
is smooth but at (1 — j)/3^ < 2 there are jumps of slopes 
in Fig. El 

Low dissipation regime and parameters M ~ 35 and 
/3 ~ 15 correspond to reality in experiments with SQUIDs 
and fit the developed theory. It is more convenient in 
experiments to obtain a set of curves as in Fig. [5] since 
usually measurements are run at a fixed temperature. 
A dependence on temperature, as in Fig. [SJ also can be 
obtained. This would provide an experimental check of 
the predicted dependences on temperature and current 



VIII. CONCLUSION 

Quantum tunneling in an asymmetric (with strongly 
different capacitances) SQUID is studied. Since capaci- 
tances play a role of masses one phase, related to a large 
mass, becomes "heavy" and remains always a constant in 
a tunneling process. Tunneling in an asymmetric SQUID 
becomes one-dimensional with a condition of optimiza- 
tion of tunneling probability with respect to a value of 
the "heavy" phase. An unusual temperature dependence 
of the tunneling probability is obtained. It has a finite 
slope at zero temperature and a transition between ther- 
mally assisted tunneling and pure activation can be not 
smooth depending on current through a SQUID. 
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